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, Abstract. The kinetic motion of the stars of a galaxy is considered within the 

framework of a relativistic scalar theory of gravitation. This model, even though 
unphysical, may represent a good laboratory where to study in a rigorous, mathematical 
^ , way those problems, like the influence of the gravitational radiation on the dynamics, 

which are still beyond our present understanding of the physical model represented by 
■ the Einstein- Vlasov system. The present paper is devoted to derive the equations of 

I the model and to prove the existence of spherically symmetric equilibria with finite 

' radius. 
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1. Introduction 

The non-relativistic kinetic model for a collisionless ensemble of particles interacting 

through the gravitational forces that they generate collectively is governed by the 
Vlasov-Poisson system: 

dtf + v-VJ-V,U-VJ = 0, (1.1) 

AU = 47rp, (1.2) 

p{t,x) = / dvf{t,x,v). (1.3) 



The Vlasov-Poisson system can be used to describe the motion of the stars in a 
galaxy, thought of as point particles, when the relativistic effects are negligible. In 
this application, the function f{t,x,v) in the previous equations is non-negative and 
gives the density on phase space of the stars in the galaxy, where t G M denotes the time 
and X, f G are the position and the velocity respectively. The function U = U{t, x) 
is the mean Newtonian potential generated by the stars. 

When the relativistic effects are not negligible, the motion of the stars has to be 
described by the Einstein- Vlasov system. In this case, the dynamics of the matter is still 
described by a non-negative scalar function /, but the unknown of the field equations 
is now a symmetric tensor of second rank, the metric of the space-time. 
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The Einstein- Vlasov system is much more comphcated than Vlasov-Poisson. The 
greatest difficulties arise from the character of the Einstein equations, which are 
essentially hyperbolic and highly non-linear even in the absence of sources, from the 
equivalence of all the coordinate systems in general relativity and from the fact that in 
the Einstein theory of gravitation, the space-time is not given in advance but is itself part 
of the solution of the equations. All these features make it extremely difficult to obtain 
global solutions of the Einstein equations coupled to any kind of matter. For instance, 
while global solutions of Vlasov-Poisson are known to be launched by general data of 
the Cauchy problem (see [HI), in the case of Einstein- Vlasov this result is known only 
for spherically symmetric solutions with small initial data (see [Zj). In particular, the 
assumption of spherical symmetry rules out one of the most interesting new features of 
general relativity, which is the propagation of gravitational waves, since this restriction 
forces the space-time to be static outside the support of the matter. It is thus not 
surprising that a rigorous mathematical theory of the gravitational radiation in general 
relativity is at present still missing, despite of the fact that its detection is currently one 
of the most challenging goals of experimental physics. 

In this paper we propose a simple, unphysical model in which the dynamics of 
the matter is still described by the Vlasov equation but where the gravitational forces 
between the particles are now supposed to be mediated by a scalar field. In spite of 
the very simplified assumptions, this model is not by any means trivial, even when it 
is further reduced to the spherically symmetric case. In fact a scalar gravitation theory 
allows also for the propagation of spherical gravitational waves. (This fact turns out to 
be very useful also for numerical investigations of the radiation in a scalar gravitation 
theory, see P). For this reason, this model may represent a comparatively easier 
framework where to study the effects of the gravitational radiation on the dynamics of 
a many particles system. It is important to underline however that the results obtained 
for this scalar model might not be representative in general for the physical model 
represented by the Einstein- Vlasov system. 

We will refer to this scalar model as the Nordstrom- Vlasov system, since the scalar 
gravitation theory which we use to describe the interaction among the particles reduces, 
in a particular frame, to the one introduced by Nordstrom in [i] (see remark 1 below). 
The equations will be derived in the next section by appealing to the more general case 
considered in |2] . In section 3 we give a few preliminary results on the Nordstrom-Vlasov 
system for which the proof traces back to the case of the Vlasov-Poisson system. All 
these results concern the spherically symmetric stationary solutions and are obtained 
by adapting the ideas of [H |Hj. 

For notational convenience, we study the model only for a single species of particle 
with unit mass (the generalization to the case of a mixture is straightforward). We also 
set the gravitational constant and the speed of light equal to one and remove the constant 
4:71 from the equations. The convention for the signature of the metric is — h -|--|-. This 
affects the definition in Cartesian coordinates of the wave operator in Minkowski space. 
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which becomes 

□ = -d^ + A. (1.4) 
2. The Nordstrom— Vlasov system 

The model which we are going to present describes the coUective motion of colhsionless 
particles interacting by means of their own self-generated gravitational forces, under the 
condition that the dynamics of the gravitational field is described in accordance to a 
simple scalar gravitation metric theory 

By scalar gravitation metric theory we mean a theory in which the gravitational 
forces are mediated by a scalar field and the effect of such forces is to induce a 
curvature in the space-time. Moreover it is assumed that the scalar field modifies the 
otherwise fiat metric only by a rescaling. Therefore the physical metric in this theory 
will be conformally fiat, that is given by 

g^, ^ A\<p)r)^,, (2.1) 

where A is a strictly positive function and, adopting Cartesian coordinates, rif^,, = 
diag(— 1, 1, 1, 1) (the use of a Cartesian frame is clearly a natural choice, due to the 
presence of a background flat metric). The metric g^i^ is assumed to be the physical one 
in the sense that it is the metric to which the matter will be coupled. The condition 
that the particles make up a colhsionless ensemble in the space-time is carried out by 
requiring that the particle density be a solution of the Vlasov equation 

dtf + ^d^^f - Tlf^d,af = 0, (2.2) 

where F^^ are the Christoffel symbols of the physical metric, / is the particle density 
and p° is determined by p° (a = 1, 2, 3) according to the relation g^uP^p''' = — 1, which 
expresses the condition that the four momentum p^ lies on the mass shell of the physical 
metric (the Creek indexes run from to 3). The previous equation is equivalent to 
postulate that the function / is constant on the geodesic flow of the metric (2.1). The 
symbols F^^ for the metric (2.1) arc given by 

where a{(f)) — ^4"^^. Hence, the Vlasov equation in a conformally flat space-time takes 
the form 

dtf + ^9,a/ - am2{p'^d,<t>)^ + ^^i^cf<p]d,af = 0. (2.3) 

We remember at this point that two tensorial quantities play an important role for 
the Vlasov equation in a curved space-time. They are the current density, 

N^ = - [ —VdP'f (2-4) 
J Po 

and the stress-energy tensor, 

T^"^ - f ^VQP^P^'f. (2.5) 
J Po 
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Q denotes the determinant of the metric and therefore, for the metric (2.1), we have 
y/g = v4^(0). The importance of A^'^ and T^'^ is due to the fact that they are both 
divergence free: 

V^T'^" = 0, (2.6) 
V^iV^ = 0, (2.7) 

where is the covariant derivative of the metric Quiy. The relations (2.6)-(2.7) hold 
independently from the dynamics of the physical metric. They are consequences of (2.2) 
alone (see [3]). 

In order to define completely the dynamics, we need now to specify the coupling 
function A{(j)) and to postulate an equation for (p. To this purpose, we remember that 
in a second order theory where the gravitational effects are mediated both by a tensor 
field and a scalar field, the field equations can be put in the form (see j2], equations 
(2.8)-(2.9)): 

i^r = 2Tr + 27(0) (9>9> + ^9^9^8,(1)0^(1)) - 25(0)^r, (2.8) 
□,.0 + j{4>)9!:''d,(Pd,^ = -a{4>)T, + (2.9) 

where E"*^ is the Einstein tensor, B the self coupling potential of the field 0, 7 an 
arbitrary function of and D^* the wave operator with respect to the metric 9*^^ . All 
the quantities in the previous equations with the suffix * are expressed in terms of the 
metric 9^^, also known as the Einstein metric. However, the Einstein frame is not the 
physical one. The physical metric, or Fierz metric, is related to 5'*^, by the conformal 
transformation 9^1, = y4^(0)(7*^,, whereas the stress-energy tensor in the Fierz frame is 
given by T^'^ = A^^{(j))T^'^ . The former rescaling law follows by the definition of the 
stress energy tensor (independent of the frame) as the variation with respect to the 
metric of the matter action (see equation (2.16) of 2J). The advantage of using the 
Einstein frame instead of the Fierz frame is due to the fact that the dynamics of the 
metric and of the field are separate from each other in the former, whereas in the latter 
they turn out to be entangled. 

Equations (2.8)-(2.9) actually define a general class of Tensor Scalar theories, since 
they contain the arbitrary function 7(0). Among these theories, a special one is singled 
out by requiring the property of scale invariance. This means that there exists in the 
theory a continuous one parameter symmetry group whose action consists in rescaling 
the coupling function A(0) by a constant factor. This assumption determines completely 
the arbitrary functions of the theory (except a constant which we fix equal to one) 
according to the relations 

A(0)=e^ 5(0) = 7(0) =0 

(see 12], equations (2.19a-b-c)). In this special theory the equation for the scalar field 
takes the form 

Dg,(j) = -T,. (2.10) 



Spherically symmetric steady states of galactic dynamics in scalar gravity 



5 



Let us explain now what we mean by simple scalar gravitation theory: we assume 
that the scalar theory satisfies the scale invariance property. Motivated by the more 
general case indicated above, we then replace in equation (2.1) with e'^. Therefore 
the Fierz metric in this scalar theory is given by 

g^^u = e'V (2-11) 
and y/Q = e^*^. Accordingly, equation (2.3) becomes 

dtf + ^9../ - Wd,<P)^ + —^dy]d,af = (2.12) 
and the mass shell condition reads 

/ = Ve-'^ + ^abPV- (2.13) 
We now introduce, beside the tensor (2.5), the "unphysical" stress energy tensor 

rp^^ = g6<^2-M^ ^ T, = e^'^T (2.14) 

and postulate for the equation 

□0 = -r„ (2.15) 

where □ is the wave operator in the fiat space-time, given by (1.4). 

Remark 1 In |4jj, the Finnish physicist Gunnar Nordstrom proposed a relativistic 
scalar theory of gravitation which was based on the field equation (2.15). For this 
reason from now on we will refer to our model as the Nordstrom- Vlasov system (or NV 
system for short). 

In the Fierz frame and for a stress-energy tensor given by (2.5), equation (2.15) 
reads explicitly 

□0 = 6^^1 ^f{t,x,p). (2.16) 

The Nordstrom- Vlasov system corresponds to the set of equations (2.12), (2.13), 
(2.16). Even though this system describes perfectly the model that we have in mind, it 
is convenient to rewrite it completely in the Einstein frame. We notice indeed that the 
particle density is still defined on the mass shell of Qf^^- In order to remove even this 
last connection with the Fierz frame, we rescale the momentum as = e'^p'^ and define 
the "unphysical" particle density by 

/.(t,x,p.) = /(t,a;,e-V)- (2.17) 
The Nordstrom- Vlasov system in the Einstein frame then becomes 

0<P = e'^l ^f4t,x,p.), (2.18) 

pI = yrr^^^, (NV,) (2.19) 

dtf. + - \iP'*d,(Pp: + d-(P)dpaJ, = 0. (2.20) 

P* P* 
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To conclude this section, we derive from (2.6) and (2.7) the conservation of energy 
and of the total number of particles. Unlike Einstein's general relativity theory in fact, 
the scalar gravitation theory under discussion here offers a natural definition for the 
local energy, in virtue of the presence of the background metric r;^^. By expanding 
equation (2.7) we get 

e-^^a^e^^AT'^ = 0. 

The previous relation shows that the quantity dt{e^'^N'^) is a pure divergence and 
therefore 

dt j dxN'^e^'^ = 0. 

Using (2.4) we then have 

jdxj dp e^'^/(t, X, p) = constant. (2.21) 

The identity (2.21) expresses the conservation of the total number of particles. The 
same relation in the NV* formulation of the model takes the form 

N^^ J J dp*e^'^/*(i,x,p*) = constant. (2.22) 

Likewise expanding (2.6) we get 

d^T^" + 6df,(/)T^"' - e-^^d^cpT = 0. (2.23) 
By means of (2.14) and (2.15), 

T = e-^^n = -e-^'^n0. 

Substituting into (2.23) and using the identity 



a>n0 = d^{d''(j)d''(j) - l-r]^"'df(j)dp(j)), 



we get 



e-^'^9^(e^'^r'^'' + dyd''(f) - \r]>'''d''(t>dp(f) = 0. (2.24) 
For i/ = 0, the previous equation leads to the conservation law 

which, using (2.5), imphes 

g{t) ^ J J dpe^V/ + dx[{dt(j>f + (V0)^] = constant. (2.25) 
For the NV* system the conservation of energy reads 

£4t) = J ^P^^^'^pIU + \j MiMf + (V0)'] = constant. (2.26) 
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3. Stationeiry spherically symmetric models 

We begin our investigation of the NV model by looking at the stationary spherically 
symmetric solutions. For this purpose we will work with the NV* formulation. For 
notational convenience we rewrite below the equations removing the suffix * (in order 
to not generate confusion we continue to call this system NV*). 



+ A0 = e^^ 



, ^^-^ [ —^==f{t,x,p), (3.1) 



dtf + • dj - m<P + ^=^)P + • d,f = 0. (3.2) 

V 1 +P V 1 + P V 1 + P 

(In this section we denote by p the vector p — (pi,P2,P3) and put = |pP)- 

Spherical symmetry means that = (f>it,r) and f{t,x,p) — f{t, Ax, Ap) for every 

orthogonal 3x3 matrix A and x,p e R^. With this restriction, NV* becomes 

-^+^^('^'0=^'V(t,r), (3.3) 

dtf + -yf= ■ dJ - [(0 + ^^=^)P + -^^] • d,f = 0, (3.4) 
f dp 

/^(^'0= / J {t,x,p), (3.5) 

I + p^ 



where "dots" and "primes" denote the temporal and the radial derivatives respectively. 
We denote by NV*[r] the system of equations which defines the stationary (i.e. 
independent of t) solutions of (3.3)-(3.5): 

J_^(0V2) = e^V(r), (3.6) 

p-dj- ^[(x ■ p)p + x] ■ dpf = 0, (NV*[r]) (3.7) 

r 

^^ir)= I -^=f{x,p). (3.8) 

For later convenience we recall here the definition in the Einstein frame of some 
important quantities: 

p(r) = e^'^*-''-* y dp-\/l +p'^f{x,p) (mass-energy density), (3.9) 
T^^r) = e^'^^'^ j -j==(^-Pj^ f{x,p) (radial pressure), (3.10) 

Vt{t) — -e^'^^'^'> I —^=^= f{x,p) (tangential pressure) (3.11) 

2 J a/I+p^ r 

and note that (3.6) can be rewritten as 

J_^(0V2) = p(^) _ (p(^) + 2VT{r)). (3.12) 
We split our analysis in three different subsections. 



Spherically symmetric steady states of galactic dynamics in scalar gravity 



8 



3.1. The Jeans theorem 

In any stationary solution of a stellar dynamics model, whether it is the Vlasov-Poisson, 
the Einstein- Vlasov or the relativistic scalar model under discussion here, the particle 
density is a functional of the invariants of the associated Vlasov equation. The statement 
according to which in the spherically symmetric case any of these invariants is in its 
turn a functional of the local energy and the angular momentum, so that in the end the 
particle density depends upon these two integrals of the characteristic system only, is 
known as the "Jeans theorem". It was established for the Vlasov-Poisson system (see 
PP) while for the Vlasov-Einstein model it was proved to be false (cf. jn])- (It remains 
open the possibility that a modified formulation of the theorem could lead to a true 
statement). In the case of the NV* system the Jeans theorem is easily seen to be valid. 
For we introduce the new coordinates 

r=\x\, w = e^^, F = e^'f'(\x\V -(x-p?), 
\x\ 

which have value in the set G = {{r,w,F) : r > 0, w G M, F > 0}. (Strictly speaking, 
the variable F could take also value zero. However, since F is an integral of motion and 
the invariant manifold defined by F = has zero measure, we can restrict ourselves to 
study the motion on G). Since / is spherically symmetric, it can be written in function 
of r, w, F, i.e. 

/ = /(r, w,F), r > 0,F > 0,w eR. 
The equation (3.7) in these new coordinates takes the form 

wdrf + {^-e''U')d^nf = 0. (3.13) 

The characteristic system for the previous equation has the form considered in ||Tj (see 
(2.11), pag. 162) with m(r) = e^'^(0'r^). Provided that is a continuously differentiable 
solution of (3.6), the function m(r) is continuous and monotonically increasing for r > 0. 
Thus we can apply the result of [TJ and conclude that any measurable function which is 
constant on the characteristics of (3.7) has the form f{r,w,F) = ^{E,F), where 

E{r, w,F) = ^w^ + ~ + U{r) = ^eV + U{r) 

and U{r) is any primitive of the function g{r) = e^'^cj)', say f/(r) = ^e'^'^- With regard 
to NV*[r], this leads immediately to the following 

Theorem 1 Let(j) G solves (3.6) (in a weak sense) and f he measurable and constant 
on the characteristics of (3.1). Then there exists a unique $ : — > [0, +oo) such that 

f{x,p) = ^E,F), (3.14) 



where E is the local energy of the particles, i.e. E = c'^a/I +p^, and F the modulus 
squared of the local angular momentum, F = e^^\x A pp. 
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Clearly, if $ is C^, then / solves (3.7) classically. 

Substituting (3.14) into the definitions (3.9)-(3.11) we get, after a transformation 
of variables, 

p{r) = -l dEE' dF ^ ^ ' ^ , (3.15) 

'P(r) = — dE dF $(E, F) E'^ - F/r^ - e^-/-, (3.16) 

p^(r) = / / dF , ^"^y^^n . (3.17) 

A direct computation leads to the following 
Lemma 1 Assume (f),p,V,VT G C^([0, R)), R> 0. Then the following equation holds: 

V'{r) = -e^'^(t)'p{r) - ^(P(r) - Vrir)), Vr G [0, R), (3.18) 

which is know as "Tolman-Oppenheimer-Volkov equation" of a spherically symmetric 
relativistic fluid hall (cf. flT] /. equation 3.2.19). 

3.2. Existence of spherically symmetric steady states 

In this subsection we prove the existence of global solutions of NV*[r] (see theorem 2 
below). We split the proof in four steps. 

Step 1. The reduced problem 

Since for any solution of NV*[r] the particle density has the form (3.14), the problem 
of proving the existence of spherically symmetric steady states can be reduced to the 
analysis of a single ordinary differential equation. For this purpose we substitute (3.14) 
into the right hand side of (3.6). After a transformation of variables we get 

A(0V2) = i7^(r,e^), (3.19) 
dr 

where 

H^ir,u)=7Tu^ I dE I dF "^l^'^^ 



v^F2 - F/r2 - m2 

The recipe to find solutions of NV=K[r] now consists in the following steps: fixing 
the functional solving (3.19) and defining / = $(F, F). Clearly, in order to carry all 
this out, some restrictions on $ must be prescribed. 

Step 2. Restrictions on $ 

The first restriction on $ is dictated by the physical requirement of finite mass of the 
corresponding steady states. By (3.15) we have 

M = An j drp^ry"^ 

= r dr r^-"" dF r ,E^^m£^. 0.20) 



Jo Je'l> 



v/F2-F/r2-e2'^' 
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The following lemma corresponds to theorem 2.1 of [S]: 

Lemma 2 Let $ : — [0, c)o) be measurable and (p be a solution of (3.6) with 
0(0) = 00- Then M < oo implies (poo = linir^+oo < oo, ^{E,F) = a.e. for 
{E,F) G (e'^°°,+oo) X [0, +oo) and the energy is finite: 



£ = M + J drr\(t)'y < oo. (3.21) 
Proof: Integrating (3.6) we get 

roc J /"OO 

/ dr— (0V2)=/ drrVV(r)<M. 
Jo dr Jq 

Since (f)'{0) = for any solution of (3.6), this implies M > lim^^+oo Since 
the function (0'r^) is increasing, then the limit in the right hand member of the latter 
inequality exists and by the condition of finite mass has to be 

lim 0'r^ < oo. (3.22) 

Thus in particular limr^+oo 0' = and since (p is increasing this implies 0oo < oo. 
Furthermore by (3.20) we have 

M>47r2 / dEE^ / dF^E,F) / dr 

Je't"^ Jo J ^F/(E'^-c^'t>oo) V £'2 — e^"*^ 

\i E > e'^°° , the integral with respect to r in the right hand side of the previous inequality 
diverges and so the mass would be infinite. Finally, (3.22) implies 0'(r) = o(r^^^'^), as 
r — s> +00, for e G [0,1) and therefore r^(0')^ ~ r^"^^ = o(l/r), for e > 1/2 and this 
proves (3.21). □ 

Remark 2 Given a solution (/, 0) with finite mass, the new solution (/, 0) defined by 
/ = e'^'l'°°f, (j) = (f) — 0OO gives a steady state which is asymptotically flat, i.e. satisfying 
the boundary condition limj.^oo 4>{r) = 0. We notice that, unlike the Vlasov-Poisson 
and the Einstein- Vlasov system, in order to obtain an asymptotically steady state from 
a solution with finite mass it is necessary to rescale the particle density (cf. the remark 
at the end of section 2 in |H]). 

Step 3. Global existence: set up 

In this step we give the basic set up to prove, in step 4, global existence and uniqueness 
for (3.19) when $ has the form: 

^E, F) = ^{E)F\ E>0, F>0,k> -1/2, 

^ G L°°((0, oo)) is non-negative, (3.23) 
3^0 > s.t. ^(E) = for E > Eo, 

which was considered in [HI E] • We start by rewriting the relevant equations according 
to (3.23). 

With $ given by (3.23), equation (3.19) becomes 

(r^0')' = vrc,,_.r2'=+V^V.(e*), (3.24) 
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where 

/■oo 

hmiu)= dEm{E){E'^ -u^)'^, m>0, (3.25) 

J u 

ca, = f s\l - sfds = + a > -1, 6 > -1, (3.26) 
Jo T{a + b + 2) 

where F is the gamma function. Moreover, the quantities (3.15)-(3.17) take the form 

p(r) = 7rr%_i^?,^.(e^), (3.27) 

p(r) = 7rr%.V|(e^), (3-28) 



TT 



^T(r) = -r^Vi,-iVf(e"). (3-29) 
where 

POO 

gm{u)= / dEm{E)E\E^ -u^r = h^+i+u^hm. (3.30) 

J u 

We also notice that in virtue of the identity xT{x) = r(a; + 1), we have 

p^(r) = (A; + l)P(r) (3.31) 
and therefore (3.12) and (3.18) become respectively 

1(0V2)' = p(r) - {2k + 3)P(r), (3.32) 

2k 

V'ir) = -e^^^'p(r) + —Vir). (3.33) 

r 

The following lemma states the regularity properties of the functions gm, hm defined by 
(3.25), (3.30) and is proved in 0^. 

Lemma 3 If m > —1 then g„i, km G C0((0, +oo)). For m>0, g^, G ^^((0, +oo)) 
and Wu > 0: 
h'^ = -2muhm-i, g'^ = -2mugm-i- 

We notice that by lemma 3, the quantities p,V and Vt are in [0, i?) provided 
that is in [0, R). The latter is then a sufficient condition for the validity of (3.33) 
in [0,i?). 

Step 4- Global existence: proof 

Here is our main result about the existence of spherically sjTiimetric stationary solutions 
of NV,: 

Theorem 2 Let 0o = 0(0) G M and $ of the form (3.23) he given. Then NV^[r] has a 
unique global solution (0, / = $(-E, F)). If the function ^ in (3.23) is , then f solves 
the stationary Vlasov equation (3.1) in a classical sense. 

Proof: We have to prove global existence for (3.24) or, equivalently, for (3.32). The 
former will be used to prove local existence, the latter to extend this solution globally. 
Integrating once (3.24) we get 



TT 
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For the function u — e'^ we then obtain the equation 



TT 



/■r 



Integrating the previous equation with the initial condition uq = ii(0) = e''*'" > we get 



u{r) = Mo + TTC 



2 



fd.^ rdrrW(r)V.(«(r)) 

s ^0 



{Tu){r). 

Now consider the set M5 = : [0, S\ [0, +00) continuous : mq < < Wo + 1}- 
By a simple calculation one can show that, for b small enough, the operator T maps 
onto itself and that T : M5 Ms is a contraction with respect to the norm. 
This gives local existence for (3.24). Standard techniques permit now to extend to a 
maximal solution on an interval [0, i?), i? > 0. \i R < 00 then 



Assume 



Then 



lim 0(r) = 00 and/or lim (j)\r) — 00. 

1 — *R~ 1 — >R~ 



lim (f){r) = a < 00. (3.34) 

r-*R~ 

J"" dr ^i<P'r'') = TTCfc^.i dr r^'^+^e^'^^'-^/i.+ i (e"^) 

<vrc,,_ii?2'^+-^e2"/.,+ .(e^«), 
where we used the fact that hm is decreasing and is increasing. Hence, 

lim (f)' < CR^^+^ < 00, 

r^R- 

where C is a positive constant. Thus it is sufficient to prove (3.34). For this purpose, 
we notice that by (3.32) we have (0'r^)' < r^p(r). Since 0' > 0, integrating twice this 
inequality we get 



0(r)<0o+ / 4 / drrV(r) 
Jo ^ Jo 



< 00 + ^Cfc -1 

<0o + Cr2'=+2, 
by which (3.34), and then global existence, follows. □ 

Remcirk 3 Since gm{u), hjn{u) — 0, for u > Eq and is increasing, then for 0o > log -E'o, 
the unique solution of theorem 2 is given by 0(r) = 0o, Vr > and f — a.e. 
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3. 3. Steady states with finite radius 

In this section we face the problem whether, among the solutions of theorem 2, there 
are steady states with compact support, i.e. such that 

E (0, +oo) : p(r) = 0, for r > R. (3.35) 

We start by pointing out that a steady state with finite radius has also finite mass (and 
so finite total energy, cf. lemma 2). For 

"R rR 



pin. pu 
M= drr^p{r) =7ic^_i / dr r^^+^^^^i (e'^ 
Jo 'wo ^ 



< 'n-Cu ^igk,i(e'^°)— < oo. 

- k, 2^/0+2 V ^2fc + 3 

Since by (3.27), p(r) = if and only if e'^^''-* > Eq, then the problem of proving 
(3.35) is equivalent to show the existence of i? G (0, +oo) such that e'^*^'-' = Eq. For this 
purpose we will need later the following simple lemma: 

Lemma 4 Let (p E and [0,R) be the maximal interval such that e^^^"* < Eq. Then 
lim e'^(") = Eq. 

r^R- 

Proof: If R < oo this is obvious. For R = +oo, let e"^"" = hmr^oo e"'^*^''^ < Eq and assume 
Q<l>ao ^ jt;^^ Then by (3.24), (r^0')' > Cr^'^"'"^, which implies, upon integration, 0' — > cx3, 
and so (f) oo, for r —>■ +oo, a contradiction. □ 

We will prove the finite radius property for steady states of the form (3.23), where 
\l/ satisfies the additional condition 

^{E) = c{Eo - E)1 + 0((Eo - E)f '), as E ^ E, , (3.36) 

where (■)+ denotes the positive part, c,6 > and 

/x>-l, k>-^, (3.37) 



+ k + fi + k + 5/2' 

This class was considered in j^, except for the condition (3.38) which is different in 
[H]. We will prove our result following the argument of 0. The result is the following: 

Theorem 3 Let (/, 0) be a solution of NV.^[r] of the form stated above. Then this 
steady state has finite radius. 

Proof: We consider only the non trivial case 0o < log-Eo (cf. remark 3). Define 
q{r) = p{r) - (2A; + 3)P(r), 



m{r) 



/ dss^q{s), r]{r) = logEQ- 
Jo 



We notice that r]{r) > for r < R and 77 ^ for r ^ R~ , where R is defined as in 
lemma 4. We also introduce the functions: 

m(r) 2?^(^) 
x{r) = ^^, yir)=r ——. 
rri[r) P[r) 
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Since rrj' = —rjx, then 

rx' = —X + + a{r)y, (3.39) 

where 

«(r) = — = {2k + 3)'' (3.40) 

m r?e2'?^/ifc+i(e<^) ^ ^ 

Using (3.27), (3.33) and lemma 3 we also get 

ry' = {2 + 2k)y - /3{r)xy, (3.41) 



where 



Pir)= - {2k + 3)77e^'^ ) . - 2{2k + 3)r/ 



+ 2(2fc + l)r/ ^ ^ . ' . (3.42) 



We reduced the field equation to the form (3.39), (3.41), which is the same stated in 
lemma 3.2 of [H] with 71 = 72 = 0. Thus if we prove that 

A) tto = infre(o,iJ) "(?") > 

B) /5o = lim,^^- /3(r) G (0,2 + 2fc), 

we conclude that R < 00. 

By means of (3.36) we have 

^{E) = c.E{El - EY^ + EO{{El - £;')f '), as ^ E^ . 

We may assume c,,, = 1, because this constant cancels in a and (3. Using the relation 

E{El - E'riE' - u'fdE = ^CaAE'o - u'r'^\ 

which is valid for u < Eq, a > —1, 6 > — 1 we get 

hmi^'^) = \c,,me^+"^+' + ©(e'^+^+^+i), E = El- e"^ 0+. 

Moreover, since 77 = log(l + e^'^e/ {Eq + e"^)), then 
1 



V=7^^e + 0{e'), e^0+. 



Thus, since e — 0^ for r — > i? (see lemma 4), we have 

C 3 

lim a{r) = {2k + 3)'^^^ 
= (/i + A; + ^)-i>0. 

Since a > on (0, i?) and limr^o+ Q^('") > (because of our assumption <^o < ^ogE^), 
then A) is proved. To prove B) we notice that 

^^(6*^) = hm+1 + e^^hm ~ e^^h^ 

= ^C^,n^£'^+'"+' + ©(e^+'^+'+l), 6 ^ 0+. 
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This immediately implies that the second term of (3.42) tends to zero as £ — > 0+ and 
therefore we get 



which lies in the required interval (0, 2k + 2) by our assumptions on Eq, /i, k. □ 
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